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I will describe recent results from the HPQCD collaboration using a new very accurate method for charm quarks in 
lattice QCD, that we have used in calculations including the full effect of u, d and s sea quarks. Multiple values of 
the lattice spacing and of the u, d and s sea quark masses allow us to extrapolate reliably, with a full error budget, 
to the real world. This opens up the field of charm physics to precision lattice QCD tests. So far we have calculated 
00 the D and Ds meson masses to 6 MeV, having fixed the charm quark mass from the r]c meson. Our D and Ds decay 

■ constants (determined to 2%) make an interesting comparison to CLEO-c results as we await improved experimental 

^ ' errors. We are also able to determine the charm quark mass to an accuracy of 1% using charmonium correlators and 

high-order continuum QCD perturbation theory. Future calculations arc briefiy discussed. 
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O 1. INTRODUCTION 
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Lattice QCD is now established as a precision tool for 'gold-plated' hadron physics enabling us both to test 

QCD to high accuracy but also to use it to determine parameters of the Standard Model such as quark masses and 

I CKM matrix elements. 
Qh' 

^ ' Figure [T] shows the current status of HPQCD lattice calculations of the masses of gold-plated mesons. These 
j-^^- are mesons which are far from strong decay thresholds and whose masses can then be accurately determined both 
theoretically and experimentally. The lattice calculations are done including the full effect of sea u, d and s quarks 
I (using the MILC collaboration gluon configurations) and at multiple values of the lattice spacing so that systematic 

■ errors from working on a space-time lattice can be removed as far as possible, and a full error budget worked out. A 
lot of the gold-plated mesons contain valence c or 6 quarks. Whilst work continues on T, i? and related mesons there 

[ has been a lot of recent progress on mesons containing c quarks, and I report on that here. One of the key issues in 
charm physics is to test lattice QCD against experiment as a precursor to using accurate lattice QCD results for B 

■ physics along with experiment to determine key CKM elements. We now have a 2% accurate lattice QCD calculation 
OO . of the decay constants of the D and Dg mesons Q and the recent experimental updates on this from CLEO-c have 

thrown up an interesting picture. This is discussed in section 2. Further theoretical tests of our charm methods have 
led to a 1% accurate determination of the charm quark mass, using methods pioneered in continnum QCD [3^, and 
this is discussed in section 3. 
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2. CHARM QUARKS IN LATTICE QCD 



In continuum QCD both charm and bottom quarks are described as heavy because their masses are significantly 
greater than the QCD scale, Aqcd, of a few hundred MeV. Special techniques can be applied to heavy quarks because 
they are nonrelativistic in their bound states and asimg) is relatively small. This is important information that can 
be applied to handling b quarks in lattice QCD also. For c quarks it turns out not to be nearly so useful because the 
charm quark mass in lattice units, rrica, is not that large for values of the lattice spacing currently in use, and will 
get even smaller as we make finer and finer lattices. We believe that a better strategy for c quarks is to treat them in 
the same way as light u, d and s quarks. This allows us to make use of light quark symmetries of the (lattice) QCD 
Lagrangian so that, for example, the annihilation rate of pseudoscalar mesons can be directly calculated without the 
need for a renormalisation constant. We can also apply standard variance reduction technqiues and find then that 
lattice statistical errors for D/Dg mesons are as small (less than 1%) as those for tt and K. 
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The gold-plated meson spectrum from lattice QCD - HPQCD 2008 
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Figure 1: The gold-plated meson spectrum from lattice QCD. Lines represent experimental values. Open circles are lattice 
results, distinguishing those hadron masses used to fix the quark mass parameters and those which were predictions ahead of 
experiment 



The worry, however, with c quarks, is that mca is still quite large and this means that nica will set the scale for 
discretisation errors that appear as a result of having a finite lattice spacing. Although these can be extrapolated 
away with multiple values of the lattice spacing, a large extrapolation will lead to a large resulting error. Typically 
we would expect, for some quantity m, to obtain a result at non-zero lattice spacing, a, that behaves as: 

m — ma=o{l + A{mca f + B{mca)^ + ...). (1) 

For a « O.lfm, {rrica) « 0.4 and (rrica)'^ terms could lead to a 20% error. These are present in most light quark 
actions (where they do not cause a big problem). We designed an action in which these terms were not present so 
that errors of only a few percent would be made even for c quarks, and once extrapolations were made to zero lattice 
spacings any remaining errors would be very small. The action is based on the improved staggered quark action used 
for sea quarks in the MILC configurations, and is called the Highly Improved Staggered Quark action (HISQ) 



2.1. Results for D and Ds mesons 

A lattice calculation proceeds by generating sets of gluon fields that are 'typical snapshots of the vacuum'. The 
MILC collaboration now have many such sets of gluon fields at 5 different values of the lattice spacing including 
the full effect of u, d and s sea quarks using the improved staggered quark formalism [t^. We have used c and 
light valence HISQ quarks on these gluon configurations to calculate correlators for rjc, D, Ds, K and tt mesons. 
By studying the correlators as a function of lattice timein the large time limit, we can then calculate the masses 
of these mesons and the decay amplitude known as the decay constant. This is defined in pure QCD as fu where 
fniTiH 0\ipJa"t5'4'\H >, and is the probability that the quark and antiquark are in the same place to annihilate 



2 



34*'' International Conference on High Energy Physics, Philadelphia, 2008 




Belle 

Charm2007 
BaBar 

hep-ex/0607094 

CLEO-c, 0806.2112, 
ICHEP08 



HPQCD HISQ u,d,s sea 
0706.1726[hep-lat] 

FNAL/MILC u,d,s sea 
LAT08 update 

ETMC u,d sea 
LAT08 prelim. 



200 225 



250 275 



300 MeV 



Figure 2: A comparison of lattice results for the D (leftmost points) and Ds (rightmost points) decay constants with experiment. 
The experimental results are obtained from the leptonic decay rate using CKM elements Vca and Vcd from elsewhere 0, . 
CLEO-c presented an update at this conference HPQCD results are from [3| using HISQ quarks. The Fermilab/MILC 
results are an update this year ETMC results are new [l^ and do not include the s sea quarks that are present in the 
real world. Neither of these results use a formalism as accurate as HISQ and hence have larger errors. Currently there is 
agreement between lattice and experiment for fo but not for , when comparing HPQCD and CLEO-c results. 



to a, W boson. When calculated in QCD, it takes into account all the strong interactions that keep the quark and 
antiquark bound into the meson. 

The quark masses need to be fixed in any QCD calculation and we do this by adjusting them until the rjc, K and tt 
meson masses are correct. The last two require extrapolation in the u/ d quark mass (the sea and valence masses are 
chosen to be the same) using chiral perturbation theory because the lattice calculation has to be done aX u/d masses 
that are larger than the physical values. D and masses have no further free parameters and we obtain values for 
them in agreement with experiment and with 6 MeV errors (again after extrapolation in the u/d mass to the physical 
point) This is a very strong test both of our control of discretisation errors and of QCD itself because it shows 
that charmonium and charm- light systems, with very different dynamics, are simultaneously described a single QCD 
Lagrangian. 

The rate for annihilation to a W{-^ Iv) is proportional to ffjV^i^ where Vat is the appropriate CKM element. Our 
results for /k/Ztt mesons can be used, with experimental rates, to determine Vus to 0.6% 2]. Our results for D 
and Ds decay constants have 2% errors and can be compared to experimental determinations obtained by dividing 
the leptonic rate by known values of Vcs and Vcd- Following our calculations the results from the experiment have 
become very exciting this year, as shown in Figure [2] A new result for fu from CLEO-c 10[ agrees well with our 
result, but values for f^^ do not. The discrepancy there amounts to 3a where a comes from experiment because our 
error is so small. This has led to speculation of new physics 1J|, since the disagreement is a serious one and is the 
only such disagreement from the 15 or so quantities that have now been accurately calculated in lattice QCD and 
compared to experiment, see Figure [TJ Final results from CLEO-c on fu^ should reduce the error bars further and 
either confirm, or not, this hypothesis. 

Meanwhile, we have continued to test our control of lattice systematic errors by calculating vector (e.g. J/ip) 
electromagnetic decay rates and hyperfine splittings such as J/tp — rjc and D* — Dg. These are in good agreement 
with experiment and we are finalising error budgets for them. Future work aims at few percent errors for semileptonic 
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form factors for D ^ K^ir for comparison with experiment. 



3. THE CHARM QUARK MASS 

We have developed, in collaboration with continuum QCD theorists, a new method for determining heavy quark 
masses. So far we have applied this to determine rric to 1% and the agreement between our result for the mass in the 
MS scheme at 3 GeV, 0.986(10) GeV and the result obtained by analogous continuum techniques, 0.986(13) 

GeV is another strong test of our control of systematic errors in handling charm quarks in lattice QCD. 

Lattice QCD calculations encompass energy scales from low-energies and long-distances (which we use to extract 
hadron masses) to relatively high-energies set by the inverse of the lattice spacing. These high-energies are well 
into the perturbative regime and we make use of this, for example, in the accurate determination of as [l7| . Our 
charmonium (rjc or J/V') correlators are also perturbative at short lattice times, before the ground state meson 
dominates (and we extract the ground state meson mass). We can access this perturbative region by taking 'time- 
moments', Gn — '^ti^lo)^G{t), of the correlator, G(t). Extrapolation of these to the continuum allows them to 
be compared to continuum QCD perturbation theory which should be accurate for small enough values of n and is 
known in some cases up to and including terms. The analogous continuum calculation [l^ effectively determines 
correlator moments in the vector case from R{e'^e~ — > hadrons) in the charm region. Our calculation in the vector 
{J li^) case agrees with experiment but our most accurate result, quoted above, comes from the pseudoscalar [rjc) |15|. 

The determination of rric from the more conventional lattice method of converting the bare lattice mass to the 
MS scheme using lattice QCD perturbation theory, agrees well with this result The accurate results for nic 

can be leveraged into an accurate result for light quark masses by determining ratios such as mc/ms from which 
renormalisation constants between the lattice and the continuum cancel [l^. 



Acknowledgments 

I am grateful to my colleagues in the HPQCD collaboration, and to members of the Fermilab lattice and MILC 
collaborations, for many useful discussions. Work supported by the STFC and the Leverhulme Trust. 



References 



[1] C.T. H. Davies et al, Fermilab lattice/HPQCD/MILC, Phys. Rev. Lett. 92:022001 (2004) ;hep-lat/0304004|. 
[2] E. FoUana et al, HPQCD collaboration, Phys. Rev. Lett. 100:062002 (2208) arXiv:0706.1726.; hep-lat]. 
[3] J. H. Kiihn, this Proceedings, arXiv:0809.1780[hep-ph]. 

[4] A. Gray et al, HPQCD collaboration, Phys. Rev. D72:094507 (2005) |hep-lat/050 70T3|. 

[5] I. F. Allison et al, HPQCD/Fermilab Lattice collaborations, Phys. Rev. Lett94:172001 (2005) 'hep-lat/0411027| . 

[6] E. Follana et al, HPQCD collaboration, Phys, Rev. D75:054502 (2007) hep-lat/0610092 . 

[7] C. Bernard et al, MILC collaboration, PoS LAT2007:090 2007 arXiv:0710.1118[hep-lat]. 

[8] L. Widhalm, Belle collaboration. Proceedings of Charm 2007. larXiv:0710.04201 

[9] B. Aubert et al, BaBar collaboration, Phys. Rev. Lett.98: 141801 (2007) |hep-ex/0607094| . 
[10] K. Ecklund et al, CLEO-c, Phys. Rev. LettlOO: 161801 (2008) arXiv:0712Tl75^hep-ex] ; B. Eisenstein et al, 

CLEO-c, arXiv:0806.2112[hep-ex]; 
[11] L. Zhang, CLEO-c, this Proceedings, |a rXiv:0 8 10.2328[ hep-ex] . 
[12] P. Mackenzie et al, Fermilab Lattice, MILC collaborations. Proceedings of LAT08. 
[13] B. Blossier et al, ETMC collaboration. Proceedings of LAT08, arXiv:0810.3145[ 

[14] See, for example, Dobrescu and Kronfeld, Phys. Rev. Lett. 100:241802 f2208l [arXrv:0803.051 2[hep-ph]. 



4 



34* International Conference on High Energy Physics, Philadelphia, 2008 

[15] I. Allison et al, HPQCD collaboration + K. G. Chetyrkin, J. H. Kiihn, M. Steinhauser and C. Sturm, 

arXiv:0805.2999[hep-lat], Phys Rev. D (in press). 
[16] J. H. Kiihn, M. Steinhauser, C. Sturm, Nucl. Phys. B778:192 (2007) [hep-ph/0702T03l . 
[17] C. T. H. Davies et al, HPQCD collaboration, arXiv:0807.1687>ep-lat]. 
[18] I. F. AUison et al, HPQCD collaboration. Proceedings of LAT08. larXiv:0810.0285[ hep-lat]. 
[19] C. T. H. Davies, G. P. Lepage et al, HPQCD collaboration. Proceedings of LAT08. [arXiv:0810.3548[ hep-lat]. 



5 



